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Abstract. A simulation model based on satellite observa- 
tions of monthly vegetation cover from the Moderate Res- 
olution Imaging Spectroradiometer (MODIS) was used to 
estimate monthly carbon fluxes in terrestrial ecosystems 
of Brazilian Amazon and Cerrado regions over the period 
2000-2004. Net ecosystem production (NEP) flux for at- 
mospheric CO 2 in the region for these years was estimated. 
Consistently high carbon sink fluxes in terrestrial ecosystems 
on a yearly basis were found in the western portions of the 
states of Acre and Rondonia and the northern portions of the 
state of Para. These areas were not significantly impacted 
by the 2002-2003 El Nino event in terms of net annual car- 
bon gains. Areas of the region that show periodically high 
carbon source fluxes from terrestrial ecosystems to the at- 
mosphere on yearly basis were found throughout the state 
of Maranhao and the southern portions of the state of Ama- 
zonas. As demonstrated though tower site comparisons, NEP 
modeled with monthly MODIS Enhanced Vegetation Index 
(EVI) inputs closely resembles the measured seasonal car- 
bon fluxes at the LBA Tapajos tower site. Modeling results 
suggest that the capacity for use of MODIS Enhanced Veg- 
etation Index (EVI) data to predict seasonal uptake rates of 
CO 2 in Amazon forests and Cerrado woodlands is strong. 
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1 Introduction 

Carbon dioxide (CO 2 ) is a major contributor to the planetary 
greenhouse effect and potential climate change. Altered cy- 
cles of carbon, water, energy, and nutrients resulting from the 
changes in Amazonian vegetation cover are expected to have 
climatic and environmental consequences at local, regional, 
and global scales. The Large Scale Biosphere- Atmosphere 
Experiment in Amazonia (LBA) is an international research 
initiative led by Brazil whose main science objectives include 
developing methods to quantify, understand, and model the 
processes controlling carbon cycling in the Amazon region 
and in the contiguous Cerrado region. 

Several previous studies of ecosystem modeling, namely 
Kindermann et al. (1996), Potter et al. (1998), Tian et 
al. (1998), Prentice and Lloyd (1998), Asner et al (2000), 
Houghton et al. (2000), Potter et al. (2001), and Foley et 
al. (2002), have examined how variations in climate affect 
the carbon balance of the Amazon basin. Most of these 
models suggested that the net annual flux of carbon by the 
basin is significantly correlated to El Nino Southern Oscilla- 
tion (ENSO) events. The Amazon basin was predicted to be 
a significant carbon sink during La Nina events, and a car- 
bon source during El Nino events. Moreover, these model- 
ing studies published before most data from the LBA project 
was generated concluded that major variations in the regional 
carbon balance are related chiefly to changes in precipitation. 
Generally, the average El Nino event has been drier than nor- 
mal, and the average La Nina period wetter in northern Ama- 
zonia. In southern Amazonia, both recent El Nino and La 
Nina periods have been drier than neutral conditions. 
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The launch of NASA’s Terra satellite platform in Decem- 
ber 1999 with the Moderate Resolution Imaging Spectrora- 
diometer (MODIS) instrument on-board initiated a new era 
in remote sensing of the Earth system with promising impli- 
cations for carbon cycle research. Direct input of satellite 
vegetation index “greenness” data from the MODIS sensor 
into ecosystem simulation models is now used to estimate 
spatial variability in monthly net primary production (NPP), 
biomass accumulation, and litter fall inputs to soil carbon 
pools. Global NPP of vegetation can be predicted using the 
relationship between leaf reflectance properties and the ab- 
sorption of photosynthetically active radiation (PAR), assum- 
ing that net conversion efficiencies of PAR to plant carbon 
can be approximated for different ecosystems or are nearly 
constant across all ecosystems (Running and Nemani, 1988; 
Goetz and Prince, 1998). 

Operational MODIS algorithms generate the Enhanced 
Vegetation Index (EVI) (Huete et al., 2002) as global image 
coverages from 2000-present. EVI represents an optimized 
vegetation index, whereby the vegetation index isolines in 
red and near infra-red spectral bands are designed to approxi- 
mate vegetation biophysical isolines derived from canopy ra- 
diative transfer theory and/or measured biophysical-optical 
relationships. EVI was developed to optimize the green- 
ness signal, or area-averaged canopy photosynthetic capac- 
ity, with improved sensitivity in high biomass regions. The 
EVI has been found useful in estimating absorbed PAR re- 
lated to chlorophyll contents in vegetated canopies (Zhang et 
al., 2005), and has been shown to be highly correlated with 
processes that depend on absorbed light, such as gross pri- 
mary productivity (GPP) (Xiao et al., 2004; Rahman, 2005; 
Huete et al., 2006). 

In this study, we present the results of CASA (Carnegie 
Ames Stanford Approach) model predictions of terrestrial 
ecosystem fluxes using 2000-2004 MODIS EVI inputs at 8- 
km spatial resolution. This is the first LBA regional model- 
ing study to infer variability in region- wide carbon fluxes for 
the Brazilian Amazon and Cerrado (savanna) ecosystems at 
such a detailed spatial resolution from NASA satellite data. 
Our CASA model (Potter et al., 1993, 1999, and 2003) has 
been designed to uniquely estimate monthly patterns in car- 
bon fixation, plant biomass increments, nutrient allocation, 
litter fall, soil carbon, CO 2 exchange, and soil nutrient min- 
eralization. We also validate the model for the first time us- 
ing LBA tower flux estimates for net ecosystem exchange of 
CO 2 in Amazon forests. 

2 Modeling methods and data sets 

As documented in Potter (1999), the monthly NPP flux, de- 
fined as net fixation of CO 2 by vegetation, is computed in 
NAS A-CAS A on the basis of light-use efficiency (Monteith, 
1972). Monthly production of plant biomass is estimated as a 
product of time-varying surface solar irradiance, S r , and EVI 


from the MODIS satellite, times a constant light utilization 
efficiency term (e mSLX ) that is modified by time- varying stress 
scalar terms for temperature ( T ) and moisture (W) effects 
(Eq. 1). 

NPP = S r EVI e max TW (1) 

The e mSLX term is set uniformly at 0.39 g C MJ -1 PAR, a value 
that derives from calibration of predicted annual NPP to pre- 
vious field estimates (Potter et al., 1993). This model cali- 
bration has been validated globally by comparing predicted 
annual NPP to more than 1900 field measurements of NPP 
(Zheng et al., 2003; Potter et al., 2007). The model uses 
the same e mSLX value for all vegetation types in the Amazon, 
but allows predicted light use efficiency to be regulated by 
monthly climate variations that vary across the region. Fu- 
ture fertilization effects from atmospheric CO 2 concentra- 
tions were considered to be inconsequential in our model, 
because the NPP algorithms in CASA are calibrated to cur- 
rent global estimates and we do not run the model into future 
years with elevated CO 2 fluxes. 

Interannual NPP fluxes from the CASA model have been 
reported (Behrenfeld et al., 2001; Potter et al., 2001) and val- 
idated against multi-year estimates of NPP from field stations 
and tree rings (Malmstrom et al., 1997). Our NAS A-CAS A 
model has been validated against field-based measurements 
of NEP fluxes and carbon pool sizes at multiple boreal for- 
est sites (Potter et al., 2001; Amthor et al., 2001; Hicke et 
al., 2002) and against atmospheric inverse model estimates 
of global NEP (Potter et al., 2003). The CASA model was 
shown to accurately simulate evapotranspiration fluxes in the 
Amazon region (Potter et al., 1998,2001). 

The T stress scalar is computed with reference to deriva- 
tion of optimal temperatures (T opt ) for plant production. The 
r opt setting will vary by latitude and longitude, ranging from 
near 0°C in the Arctic to the middle thirties °C, in low lati- 
tude deserts. The W stress scalar is estimated from monthly 
water deficits, based on a comparison of moisture supply 
(precipitation and stored soil water) to potential evapotran- 
spiration (PET) demand using the method of Priestly and 
Taylor (1972). Latent heat fluxes were predicted as esti- 
mated evapotranspiration (EET) amounts, which was deter- 
mined the lower value in a comparison of monthly PET and 
precipitation plus stored soil water volumes. 

The Moderate Resolution Imaging Spectroradiometer 
(MODIS) 1-km land cover map (Friedl et al., 2002) aggre- 
gated to 8-km pixel resolution was used to specify the pre- 
dominant land cover class for the W term in each pixel as ei- 
ther forest, savanna ( Cerrado ) crop, pasture, or other classes 
(i.e., water or urban area). Monthly mean surface air tem- 
perature and precipitation grids for model simulations over 
the years 2000-2004 came from NCEP reanalysis products 
(Kistler et al., 2001). Monthly mean inputs of solar radiation 
flux to the model were derived from top of the atmosphere 
shortwave radiation budget products of Laszlo et al. (1997 
and 2006). 
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Evapotranspiration is connected to water content in the 
soil profile layers (Fig. 1), as estimated using the NASA- 
CASA algorithms as described by Potter (1999). The soil 
model design includes three-layer (Mi— M 3 ) heat and mois- 
ture content computations: surface organic matter, topsoil 
(0.3 m), and subsoil to rooting depth (1 to 10 m). These 
layers can differ in soil texture, moisture holding capacity, 
and carbon-nitrogen dynamics. Water balance in the soil 
is modeled as the difference between precipitation or vol- 
umetric percolation inputs, monthly estimates of PET, and 
the drainage output for each layer. Inputs from rainfall can 
recharge the soil layers to field capacity. Excess water per- 
colates through to lower layers and may eventually leave the 
system as seepage and runoff. 

Based on plant production as the primary carbon and ni- 
trogen cycling source, the NASA-CASA model is designed 
to couple daily and seasonal patterns in soil nutrient miner- 
alization and soil heterotropic respiration (Rh) of CO 2 from 
soils worldwide. Net ecosystem production (NEP) can be 
computed as NPP minus Rh fluxes, excluding the effects of 
small-scale fires and other localized disturbances or vegeta- 
tion regrowth patterns on carbon fluxes. The NASA-CASA 
soil model uses a set of compartmentalized difference equa- 
tions with a structure comparable to the CENTURY ecosys- 
tem model (Parton et al., 1992). First-order decay equations 
simulate exchanges of decomposing plant residue (metabolic 
and structural fractions) at the soil surface. The model also 
simulates surface soil organic matter (SOM) fractions that 
presumably vary in age and chemical composition. Turnover 
of active (microbial biomass and labile substrates), slow 
(chemically protected), and passive (physically protected) 
fractions of the SOM are represented. Along with moisture 
availability and litter quality, the predicted soil temperature 
in the M\ layer controls SOM decomposition. 

The NASA-CASA soil carbon pools were initialized to 
represent storage and flux conditions in near steady state (i.e., 
an annual NEP flux less than 0.5% of annual NPP flux) with 
respect to mean land surface climate recorded for the pe- 
riod 1979-1981 (New et al., 2000). This initialization pro- 
tocol was found to be necessary to eliminate any notable dis- 
continuities in predicted NEP fluxes during the transition to 
our model simulation years of interest prior to MODIS EVI 
availability, which were run on a monthly time step starting 
January 1982 to December 2000. Initializing to near steady 
state does not, however, address the issue that some ecosys- 
tems are not in equilibrium with respect to net annual carbon 
fluxes, especially when they are recovering from past dis- 
turbances. For instance, it is openly acknowledged that the 
NASA-CASA modeling approach using 8 -km satellite data 
inputs cannot capture all the carbon sink effects of forest re- 
growth from recent wood harvest activities (Turner, 2005), 
although impacts of major wildfires are detectable (Potter et 
al., 2005). Higher resolution (250-m) MODIS EVI data sets 
are currently in the evaluation phase for use in NASA-CASA 
model runs for deforested and regrowth areas. 


Whereas previous versions of the NASA-CASA model 
(Potter et al., 1993, 1999) used a normalized difference veg- 
etation index (NDVI) to estimate FPAR, the current model 
version instead has been recalibrated to use MODIS EVI 
datasets as direct inputs to Eq.(l) above. In long-term (1982- 
2004) simulations, continuity between AVHRR and MODIS 
sensor data for inputs to NASA-CASA is an issue that must 
be addressed by recalibration of annual NPP results post 
2000. NASA-CASA model predictions with 2001 monthly 
MODIS EVI inputs have been adjusted using the same set of 
field measurements of NPP (Olson et al., 1997; Potter et al., 
2003, 2007; Zheng et al., 2003). To best match the predic- 
tions with previously measured NPP estimates at the global 
scale (R 2 =0.91), the model e max term for 2001 MODIS EVI 
inputs was reset to 0.55 g C MJ -1 PAR. 

3 Evaluation of net carbon flux results at tower sites 

Measurements of net ecosystem carbon exchange from tower 
studies of eddy correlation can be used with confidence for 
the purpose of carbon model evaluation, but only after inde- 
pendent confirmation research can show that a given tower 
flux time series can in fact yield consistently high quality 
measurements of NEP. Tower flux data from Tapajos Na- 
tional Forest (TNF) sites near Santarem, Para have been val- 
idated using radon as a tracer (Martens et al., 2004), which 
indicated that eddy flux carbon balance were realistic there. 

The TNF tower site is in the east-central Brazilian Ama- 
zon (2°54 / S, 54°57 / W). The region receives an average an- 
nual rainfall of about 2000 mm, varying between 600 and 
3000 mm annually, and suffers severe drought stress during 
ENSO episodes. Soils are deeply weathered Oxisols (Hap- 
lustox), with no concretions or impeding layers in the upper 
12 m. The water table is 100 m below the soil surface. 

Comparison of NASA-CASA predicted net ecosystem 
exchange (NEE, which is merely the negative transform 
of monthly NEP flux) to measurements at the TNF site 
(Hutyra et al., 2007) showed a close month-to-month match 
(R 2 =0.57) between the two flux estimates over the three- 
year period, 2002-2004 (Fig. 2). This comparison assumed 
a one-month lag correction in the NASA-CASA predicted 
NEE flux, due mainly to the rapid response of heterotropic 
respiration (Rh) flux of carbon to changes in soil water that 
the monthly time step in NASA-CASA cannot update on a 
shorter time interval. Nonetheless, the NASA-CASA model 
predictions closely track the seasonality and magnitude of 
tower measurements of NEE, with net source fluxes (posi- 
tive NEE) during the wet season periods and rapid transition 
to net sink fluxes (negative NEE) during the dry season peri- 
ods of each year. By separating the NASA-CASA predicted 
NEE flux into its component (NPP and Rh) fluxes, (Fig. 3), 
the opposite seasonality of forest NPP and soil Rh fluxes of 
CO 2 can be easily discerned. Predicted soil Rh fluxes were 
nearly as variable as predicted NPP fluxes over the course of 
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Fig. 1. Schematic representation of components in the NASA-CASA model. The soil profile component (a) is layered with depth into a 
surface ponded layer (Mo), a surface organic layer (Mi), a surface organic-mineral layer (M2), and a subsurface mineral layer (M3), showing 
typical levels of soil water content (shaded) in three general vegetation types (DeFries et al., 1995). The production and decomposition 
component (b) shows separate pools for carbon cycling among pools of leaf litter, root litter, woody detritus, microbes, and soil organic 
matter, with dependence on litter quality (q). 
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Fig. 2. Comparison of the 8-km resolution NASA-CASA predicted 
NEE versus measured monthly NEE from the Tapajos National For- 
est (TNF) tower flux site near Santarem, Para. 

an entire year. Rh fluxes were highest in the wet season and 
declined rapidly at the same time that dry season NPP fluxes 
increased rapidly. 

Other than the TNF tower flux data sets, we are not aware 
of similarly comprehensive field validation work completed 
for other Amazon forest tower sites. Furthermore, measure- 
ments at the ZF2 Manaus tower site (which used a multiple 
regression analysis to identify times when nighttime mea- 
surements are reliable) suggested that only rarely could the 
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Fig. 3. Monthly patterns of the NASA-CASA predicted NPP and 
Rh fluxes at the Tapajos National Forest (TNF) tower flux site. 


eddy flux carbon balance be used with confidence for daily 
carbon balance purposes. This is due to low atmospheric tur- 
bulence. Chambers et al. (2004) were able to identify only 
about lOOh out of over 4000 measurements of the nighttime 
eddy flux collected at this tower site to develop a net carbon 
budget, which would not be sufficient for long-term seasonal 
evaluation tests of ecosystem carbon models such as NASA- 
CASA. 


Biogeosciences, 6, 937-945, 2009 


www.biogeosciences.net/6/937/2009/ 


C. Potter et al.: Carbon sinks in the Brazilian Amazon 

Nevertheless, Chambers et al. (2004) directly measured 
respiration rates from live leaf, live wood, and forest 
soil surfaces to derive an indirect ecosystem NPP flux 
estimate of 900 g Cm -2 yr -1 , and further partitioned be- 
lowground respiration fluxes for an annual Rh estimation 
of 850gCm -2 yr -1 (and an NEP sink flux estimate of 
+50 g C m -2 yr -1 ) at the ZF2 Manaus tower site for the pe- 
riod 1999-2000. In comparison, our NASA-CASA predic- 
tions for annual NPP flux at this site ranged from 839 to 
867gCm -2 yr -1 for the period 2001-2004 with a mean 
annual NPP flux of 857gCm _2 yr _1 , and mean annual 
Rh flux of 846 g Cm -2 yr -1 , for a NEP flux estimate of 
+11 gCm -2 yr -1 at the ZF2 Manaus site. Taking into con- 
sideration the high levels of topographic and soil variations 
affecting forest carbon fluxes over this area of the central 
Amazon, the apparent difference of less than 1% between 
measured and our modeled annual Rh fluxes at the Manaus 
site is an important confirmation of the NASA-CASA perfor- 
mance at this additional LB A study location. 

4 Regional carbon flux results 

Annually summed NPP fluxes for the Legal Amazon re- 
gion of Brazil decreased gradually over the period of 2000- 
2004, from 4.34PgCyr -1 (1 Pg=10 15 g) in 2001 to a low of 
4.25PgCyr -1 in 2002, and then again to 4.26PgCyr -1 in 
2004. (Fig. 4a). The smoothed trend in monthly NPP rates 
over the region, which were typically predicted by the model 
to peak during the period of October to January, declined 
steadily over the period from 2001 to early 2003, and then 
recovered slowly over the subsequent two years for the re- 
gion. 

Annually summed NEP fluxes for the Legal Amazon re- 
gion varied from year to year over the period of 2000- 
2004, from a relatively low CO 2 emission source of 
— 0.07PgCyr _1 C in 2000 to a higher of CO 2 emission 
source of — 0.13PgCyr _1 in 2002 (Fig. 4b). The modeled 
NEP fluxes were controlled strongly by monthly predicted 
Rh fluxes of CO 2 from soil microbial activity, in which case 
Rh accounted for more that 32% of the variation in the pre- 
dicted NEP sink fluxes region- wide (based on the correlation 
coefficient r, significant at a level of /?<0.01). Monthly NEP 
predictions were not highly correlated (r= 0.2; not significant 
at p< 0.1) with predicted NPP fluxes on a month-to-month 
basis. 

Geographic patterns in annually summed NPP and 
NEP fluxes confirm that the year 2002 stood out from the 
other years 2000-2004 with lower ecosystem productivity 
and relatively large carbon source fluxes across the region 
(Fig. 5). This year-to-year contrast can be most readily ob- 
served in the north eastern portion of the Brazilian Amazon 
region, specifically along the entire Tocantins River basin 
as far south as the city of Maraba, along the Tapajos River 
basin located south of the Amazon River main stem between 
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(b) 2000 2001 2002 2003 2004 


Fig. 4. (a) Net Primary Production (NPP) in ecosystems of the 
Brazilian Amazon region from 2000 to 2004. Thin line is the 
monthly predicted NPP; thick line is the 12-month running average 
NPP. (b) Net Ecosystem Production (NEP) of the Brazilian Ama- 
zon region from 2000 to 2004. Thin line is monthly predicted NEP; 
thick line is the 12-month running average NEP. 


longitudes 57° and 54° W, and also in the zones located 
just north and as far west as Lago Badajos from the city 
of Manaus, mainly between longitudes 64° and 60° W. Due 
in part to a 5-15% decline in annual predicted NPP, pre- 
dicted annual NEP for these specific Amazon forest zones 
typically switched from a net sink flux on the order of +20 
to +90 g C m -2 yr -1 during 2001 and 2003 to net source flux 
on the order of —30 to —75 g C m -2 yr -1 during 2002. 

According to the climate data inputs to these NASA- 
CASA runs, total precipitation across the Legal Amazon 
states of Brazil was about 3% above the regional five-year 
(2000-2004) mean in 2000 and 2001 and was about 6% be- 
low the regional five-year mean in 2002. Zeng et al. (2008) 
analyzed the regional impacts of drought in the Amazon cul- 
minating in 2005 and attributed the causes to the combination 
of 2002-2003 El Nino event and another dry period in 2005 
that was coupled to warm subtropical North Atlantic Ocean 
conditions. According to standardized precipitation index 
(SPI) values derived from satellite rainfall data (Janowiak 
and Xie, 1999), the 2002-03 El Nino event was as strong 
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Fig. 5. Predicted annual NEP fluxes from the NASA-CASA model from 2001 to 2004 over the Brazilian Amazon region. 



0 Latitude 

Fig. 6. Correlation coefficients (r values) for CASA model monthly 
NPP and Rh fluxes of CO 2 versus monthly precipitation amounts 
(PPT) or predicted EET fluxes of water over the period 2001-2004. 
Lines are based on ten randomly selected latitude zones, starting 
from north of the Equator (positive latitude values) to south (nega- 
tive latitude values) of the Equator. Absolute values of r>0.25 are 
significant at the level of p<0.05 (n = 48 months). 


as the 1997-98 El Nino event across the Amazon, and con- 
tributed to near record-low stream flows and enhanced fire 
frequency by 2006. 


The NASA-CASA model predictions for the Cerrado 
woodland zones represented in the states of Mato Grosso and 
Goias were also impacted by the 2002-2003 El Nino event 
with lower ecosystem productivity and relatively large car- 
bon source fluxes (Fig. 5). Predicted annual NEP for these 
specific Cerrado zones typically switched from a net sink 
flux on the order of +30 g C m -2 yr -1 during 2001 and 2003 
to net source flux on the order of — 25 gCm -2 yr -1 during 
2002. 

A unique contribution of the CASA model has been de- 
tailed spatial mapping of interannual variation in ecosystem 
carbon fluxes (Potter et al., 2002). To further explore poten- 
tial climate controls over Amazon carbon source (Rh) and 
sink (NPP) fluxes, zonal correlations between aggregated 
CASA CO 2 fluxes and monthly precipitation amounts (PPT) 
or predicted EET fluxes of water were computed over the 
period 2001-2004. Results in Fig. 6 revealed that interan- 
nual variations in both precipitation and latent heat fluxes (as 
represented in EET) strongly controlled the predicted car- 
bon source and sink fluxes in zones south of 7° S latitude, 
whereas latent heat fluxes alone were far more important con- 
trollers of NPP than were monthly precipitation amounts in 
zones north of 7° S latitude. On the basis of these results, we 
hypothesize that rainfall is rarely limiting to net ecosystem 
carbon in the Amazon forest zones north of about 7° S lati- 
tude, and that seasonal variation in solar irradiance is a more 
important controller of annual NPP levels in these areas. 
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5 Discussion 

Previous ecosystem modeling studies such as those of Tian 
et al. (1998) and Foley et al. (2002) have implied that ma- 
jor variations in the regional carbon balance of the Amazon 
are related chiefly to ENSO precipitation patterns. Our cur- 
rent study results, combined with previous CASA findings 
from Potter et al. (2004) are consistent with these previous 
findings, mainly for the northern portion of the Amazon re- 
gion. Our prior model analysis over the period 1982-1998 
suggested that temperature effects of ENSO cycles can also 
have notable impacts on Amazon carbon fluxes over decade- 
long time periods. Specifically, model anomalies in NPP and 
NEP over large areas of the Amazon region east of 60° W 
longitude were strongly correlated with the Southern Oscil- 
lation Index (SOI). To explain NEP anomalies, reasonably 
strong temperature correlations with SOI were even more 
widespread regionally than were strong correlations of pre- 
cipitation with SOI. 

It is worth noting that NASA-CASA has the potential 
to generate results that are very different than other mod- 
els such those described by Tian et al. (1998) and Foley et 
al. (2002), largely because NASA-CASA relies on monthly 
satellite data as an input. The MODIS and AHVRR sensors 
generate observations of the actual land vegetation dynam- 
ics, which are lacking in the other cited models, and which 
can be sensitive to variability on temperature as well as pre- 
cipitation. 

At the regional scale, the CASA model of Potter et 
al. (2004) estimated NPP for the total Amazon region (be- 
yond Brazil) to vary between 8.7 (in 1983) and 9.8 PgC 
(in 1997) yr -1 , with a seasonal anomaly range of about 
0.1 PgC month -1 . This study also inferred an increase in 
NPP for the region over the 17-year time period of 1982- 
1998. The predicted NEP flux for atmospheric CO 2 over the 
region varied between an annual source (to the atmosphere) 
of 0.17PgCyr -1 in 1983 to a sink (from the atmosphere) 
of 0.64PgCyr -1 in 1989, with a seasonal anomaly range of 
up to 0. 1 Pg C month -1 . The NPP and NEP predictions from 
our new MODIS -era study of the Brazilian Amazon states 
fit well within these previous total flux ranges, and have the 
distinct advantage of matching tower fluxes closely. Due in 
part to improved sensitivities to declines in annual predicted 
primary production detected uniquely by the MODIS sen- 
sor data, predicted annual NEP for specific Amazon forest 
zones can be mapped out as switching from a net sink flux to 
a net source during ENSO cycles, a capability never before 
available for assessing long-term ecosystem carbon balance 
across this remote region. Our newest CASA results include 
carbon fluxes for the Cerrado region as well, which deserves 
much greater study in high-resolution (1-km or finer) model 
runs that are being designed next for CASA to include many 
different crop systems from converted savanna ecosystems. 

Inclusion of deforestation sources of CO 2 emission in 
tropical zones is a crucial next step in developing a full 


carbon balance for the region. Previously, VanderWerf et 
al. (2003) used the CASA model with inputs from satel- 
lite observations of global rainfall from the NASA Tropical 
Rainfall Monitoring Mission (TRMM) to account for the ef- 
fects of fire on regional carbon stocks and ecosystem car- 
bon fluxes. Annual carbon emissions from fires in the Legal 
Amazon were estimate to range from 0.2 to 0.5PgCyr -1 . 
In comparison, Potter et al. (2002) used the CASA model 
production and biomass predictions together with Landsat- 
derived mapping of burned areas for the Legal Amazon 
states to estimate total ecosystem source to range from 0.2 
to 1.2PgCyr -1 for the period 1992-1993. The CARbon 
and Land-Use Change (CARLUC) model estimates the net 
flux caused by deforestation and forest re-growth (Hirsch et 
al., 2004). CARLUC predicted that the net flux to the atmo- 
sphere from the Legal Amazon area during the period from 
1970 to 1998 reached a maximum of 0.35 Pg C yr -1 in 1990, 
with a cumulative release of 7 Pg C. The net flux is higher 
than predicted by an earlier study (Houghton et al., 2000) by 
a total of 1 Pg C over the period 1989-1998, mainly because 
CARLUC predicts relatively high mature forest carbon stor- 
age compared with the datasets used in earlier studies. 

6 Conclusions 

The results of NASA-CASA model simulations of NPP and 
NEP across the Brazilian Amazon and into the Cerrado re- 
gions from 2001-2004 have several noteworthy implications 
for LB A assessment studies: 

- Areas of the region that show consistently high carbon 
sink fluxes in terrestrial ecosystems on a yearly basis are 
the western portions of the states of Acre and Rondonia 
and the northern portions of the state of Para. These 
areas were not significantly impacted by the 2002-2003 
El Nino event in terms of net annual carbon gains. 

- Areas of the region that show periodically high car- 
bon source fluxes from terrestrial ecosystems to the at- 
mosphere on yearly basis are throughout the state of 
Maranhao and the southern portions of the state of Ama- 
zonas. These areas were predicted as being strongly af- 
fected by the 2002-2003 El Nino event in terms of re- 
ducing net annual carbon gains. 

- Cerrado region (woodland savanna formations) showed 
behavior of NEP on yearly basis similar to the southern 
part of the Amazon but fluxes were approximately one 
third of those measured for the Amazon region. 

We note in closing that an advantage of combining ecosys- 
tem modeling with satellite observations for vegetation cover 
properties is to uniquely enhance the spatial resolution of 
source and sink patterns for CO 2 in the terrestrial biosphere. 
Using MODIS land products, carbon modelers have begun 
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to identify numerous relatively small-scale patterns through- 
out the world where terrestrial carbon fluxes have varied in 
recent years between net annual sources and sinks. Predic- 
tions of NEP for these areas of high interannual variability 
will require further uncertainty analysis of carbon model es- 
timates, with a focus on both flux algorithm mechanisms and 
potential scaling errors to the regional level. 
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